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We calculate matrix elements of the axial current for static-light mesons and baryons in lattice 
QCD with dynamical domain wall fermions. We use partially quenched heavy hadron chiral 
perturbation theory in a finite volume to extract the axial couplings gi, g2, and ^3 from the data. 
These axial couplings allow the prediction of strong decay rates and enter chiral extrapolations of 
most lattice results in the b sector. Our calculations are performed with two lattice spacings and 
with pion masses down to 227 MeV. 
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Axial couplings of heavy hadrons from domain-wall lattice QCD Stefan Meinel 

1. Introduction 

The low-energy dynamics of heavy-light mesons and baryons can be described by heavy- 
hadron chiral perturbation theory (HH;t;PT), an effective field theory for QCD that incorporates 
both chiral symmetry and heavy-quark symmetry [1]. HH;^PT is essential for controlling light- 
quark-mass extrapolations of lattice QCD data in the heavy-quark sector (see for example Ref. [2]). 

At leading order in the heavy-quark and chiral expansions, the HH;^PT Lagrangian contains 
three axial coupling constants that determine the strength of the interactions between heavy-light 
hadrons and pions: one coupling (denoted as g\) for the heavy-light mesons, and two additional 
couplings (denoted as g2, gs) for the heavy-light baryons. These axial couplings ai^e calculable 
from QCD, and their determination enables quantitative predictions for many heavy-light hadron 
properties (such as masses, decay widths, and various matrix elements) using HH;^PT. The chiral 
loop contributions that lead to the nonanalytic dependence of such properties on the light-quark 
masses are proportional to products of the relevant axial couplings. While ^i has received much 
attention in the past because of its role for B mesons, the lesser-known couplings g2 and ^3 are im- 
portant for flavor physics with heavy baryons. The bottom baryon sector provides complementary 
information to B mesons for constraining the helicity structure of possible new physics [3]. 

The calculation of gi^2,3 from the underlying theory of QCD must be done nonperturbatively, 
and hence on a lattice. The mesonic coupling gi had been studied previously in lattice QCD with 
Uf = Oornf = 2 dynamical flavors [4]. In the following, we present a complete determination of all 
three axial couplings ^1.2,3 using rif = 2+1 domain- wall lattice QCD [5,6]. Our choice of lattice 
parameters (low pion masses, large volume, two lattice spacings) and our analysis method (fits to 
the axial-current matrix elements using the correct next-to-leading-order formulae from HH;^PT 
[7]) allow us to control all sources of systematic uncertainties. 

2. Heavy-hadron chiral perturbation theory 

We begin with an introduction to HH;^PT. This theory combines the chiral expansion with 
an expansion in powers of Aqcd/iiq, where niQ is the heavy-quark mass. We work at the leading 
order in the heavy-quai^k expansion, where the spin of the light degrees of freedom (si) is conserved 
and the heavy-quark spin decouples. The lowest-lying heavy-light mesons with 5/ = 1/2 form 
multiplets with / = and / = 1 , which can be combined into a single field H' : 

h' = [-p^y5+p;y]^-^, with in=(Xo]^ (^*'')^ = (r*0- ^^-^^ 

(We consider SU{2) HH;^PT and use the notation for bottom hadrons.) Similarly, the baryons 
with si = I form multiplets with J = 1/2 and / = 3/2. These are described by Dirac and Rarita- 
Schwinger fields B'^ and B^'^, which are symmetric in the flavor indices and can be combined into 
a single field S'^: 

nr / y+ XvO\ / y*+ _Lv*0' 



S^i = \hiY,+^,)Y5B'^+B,\ With (B'J) =\ / ^^_'^ , {B*'% - , , 
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On the other hand, the si =0 baryons (7=1 /2) are antisymmetric in flavor and include only the 
Afo in the SU{2) case: 

The leading-order HH;^PT Lagrangian, describing the interactions of the fields (2.1), (2.2), and 
(2.3) with pions, is given by 



^ = ( kinetic terms ) + gi tro [Hi{£/^'y jY^y^W] 
In the terms with g\, g2, and g^,, the pion field ^ = \/L = exp(/<I>//) appears through 



(2.4) 



■^ / 

which is an axial- vector field. 

3. Axial current matrix elements 

To determine the axial couplings g, that appear- in the chiral Lagrangian (2.4) from QCD, 
one can calculate suitable hadronic observables in both HH;t;PT and lattice QCD. The expressions 
derived from HH;^PT are then fitted to the lattice data, and in these fits the axial couplings are pa- 
rameters. The simplest observables that are sensitive to gj are the zero-momentum matrix elements 
of the axial current between heavy-hadron states. In QCD, the isovector axial current is given by 
the quark current 

.a(QCD) _T" 

V =^y7M75?- (3.1) 

The corresponding hadronic current in HH;^PT can be obtained from the Lagrangian (2.4) using 
the Noether procedure. At leading order, the relevant part of the current that contributes to the 
matrix elements reads 



^«(;,PT,Lo) ^ ^^^^ H^irlJjY.r^W -ig2e^,,,slv^{rlJj{S^) 



J^V'*: 



+ V2g,\{S,UTlJjTJ>^ + T,^{rlJj{S,) 



jk 



(3.2) 



with T?^ = 2 (i^^t"(^ + ^t"^^) . One finds the following matrix elements for A^ ~'^lt~ '"^m' 

{P*''\A^r') = -2{giUe;, 

{S'''\A^\S"') = -{i/V2){g2Uv^ea^ypU''uP, 

{S'''\A^\r''') = -{g,UU^^, (3.3) 

where e^ is the polarization vector of the /"*'' meson and '^ is the Dirac spinor of the T"^" baryon. 
For the si = \ baryons, we work directly with external 5'^ states (which contain the degrees of 
freedom of both J = 1/2 and J = IjT) and the J^'s are the corresponding generalized spinors [6]. 
At leading order in the chiral expansion, the "effective axial couplings" in (3.3) are equal to the 
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Figure 1: One-loop diagrams contributing to the matrix elements of the axial current in HH;|^PT: (a) 
wavefunction renormalization diagram, (b) tadpole diagram, (c) sunset diagram [7]. 

axial couplings in the Lagrangian, (^,)eff L q = Si- ^t next-to-leading order, the matrix elements 
receive corrections from pion loops (Fig. 1) and analytic counterterms. The NLO expressions for 
(§r)eff. both in the unquenched and the partially quenched theories, and for a finite volume, can be 
found in Ref. [7]. 

4. Lattice calculation 

To calculate the matrix elements (3.3) in lattice QCD, we use the following interpolating fields 
for the heavy hadrons, 

P' = (75 ) a/3 Qaa C/3 ' P*^ = {rii)ap Qaa C/3 ' 

S'fia = £abciCY^i)pyqi,pqiyQca, T'^ = £ahc{Cy5)pyq'^pqi^Qca, (4.1) 

where q' denotes a smeared light-quark field of flavor /, and Q denotes the static heavy quark field 
(here we set v = (1,0,0,0)). Just as in HH;^PT, the interpolating field S'^a couples to both the 
7=1/2 and J = 3/2 baryon states with 5/ = 1. We use the domain- wall action [8] for the light 
quarks, and the Eichten-Hill action [9] with HYP- smeared temporal gauge links [10] for the heavy 
quarks. To optimize the signals and analyze heavy-quark discretization effects, we generated data 
for ?iHYP = 1 , 2, 3, 5, 10 levels of HYP smearing. The final results for the axial couplings are based 
on data with «hyp = 1,2,3 only. The calculations are performed with the local 4-dimensional axial 
current, given by 

A^=ZAdaa{7n75)apUap, (4.2) 

where Za is determined nonperturbatively [11]. We compute the following ratios of three-point and 
two-point functions. 



Riit, t') 
R,{t, t') 



\ Zl=,{P*'nt)A^{t')Pl{0)) 
"3 (P"(04(0)) [ 



(4.3) 



li:l.y=i{S''^{t)A^^{t')T,M){T'"{t)A^'^{t')sl{Or''^' 



3 LU(^"'nt)s',AO)){T'"{t)T,M: 



where the source and sink hadron interpolating fields are placed at a common spatial point x be- 
cause of the static heavy quark, and we write A'^(f') = T.x"^^{f' ^'^')- I^i Eq. (4.3) we also removed 
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sea) 
u/d 



am 



u/d 



a (fm) 



m 



W 



(MeV) m^r^ (MeV) values of? /a 



24^ X 64 


0.005 


0.005 


0.1119(17) 


336(5) 


24^ X 64 


0.005 


0.002 


0.1119(17) 


270(4) 


24^ X 64 


0.005 


0.001 


0.1119(17) 


245(4) 


323 X 64 


0.006 


0.006 


0.0848(17) 


352(7) 


323 X 64 


0.004 


0.004 


0.0849(12) 


295(4) 


323 X 64 


0.004 


0.002 


0.0849(12) 


227(3) 



336(5) 


4,5,6,7,8,9,10 


304(5) 


4, 5, 6, 7, 8, 9, 10 


294(5) 


4, 5, 6, 7, 8, 9, 10 


352(7) 


13 


295(4) 


6,9,12 


263(4) 


6,9,12 



Table 1: Lattice parameters, m)^^ and m)^^ denote the valence-valence and valence-sea pion masses. 
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Figure 2: Left panel: ratios /?,(f, t') for the source-sink separation t/a = 10, along with extracted values 
Ri{t) (shaded regions). Right panel: extrapolation of /?,(f) to infinite source-sink separation. All data shown 



here are for a = 0. 112 fm, am 



val) 
u.d 



■ 0.002, «HYP = 3. 



the free spinor indices which trivially originate from the static heavy-quark propagator. By insert- 
ing complete sets of states into (4.3), one can show that /?,(?, f/2) = (gi)eff + •••, where the dots 
indicate contributions from excited states that decay exponentially with t [6]. 

Our calculations use RBC/UKQCD gauge field configurations with 2-1-1 dynamical quark fla- 
vors [11]. The main parameters of the ensembles and the domain- wall propagators we computed 
on them are given in Table 1 . For the three-point functions we use pairs of light-quark propagators 
with sources at a common spatial point x and separated by t/a steps in the time direction. As can 
be seen in the table, we have data for multiple values of t/a. Numerical examples for the ratios 
(4.3) are shown in Fig. 2 (left). Equivalently to using /?,(?, t /2), we average /?,(?, t') over t' in the 
central plateau regions, and we denote these averages as /?,(?)■ We then perform fits of the form 
^KO = (,?()eff — ^( e^^'' , as shown in Fig. 2 (right). A detailed discussion can be found in Ref. [6]. 
These fits provide the effective axial couplings (g()eff(fl^,'Mjr^^ ,'m}j^ ,nHYp) for all combinations of 
the lattice spacing, the pion masses, and the heavy-quark smearing parameter «hyp- We fit the 
data for (g,-)eff(«,"4"\"4'''',nHYp) with 



?l)eff = gi 



Jeff 



(=2, 3 



1 +/,(§2, ^3, m^r\m^r\ A, L) +C^'^ r,Jvv)i2 ^^(vs) ,Jvs),2 ^^. 



[mTT+dLn^,,a^ ,(4.4) 
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Figure 3: Fits of (g,)eff using Eq. (4.4). The plot shows the fitted functions, evaluated at mjj = »z^ = mjt 
and in infinite volume, for «hyp = 3. The baryonic matrix elements (§2,3)eff develop small imaginary parts 
below the S -^ Tn decay threshold at m„ = A, and only the real parts are shown here. The dashed line 
corresponds to a = 0. 112 fm, the dotted line to a = 0.085 fm, and the solid line to the continuum limit. The 
±1(7 regions are shaded. The data points (circles: a = 0.1 12 fm, squares: a = 0.085 fm) have been shifted to 



(vv) 



,(-)^ 



infinite volume for this plot, and the partially quenched data (m^ < m)^ ) are included using open symbols 
at m-,1 = m]^^ , even though the fit functions have slightly different values for these points. 



where the functions fj are the NLO loop contributions from 56^^(412) HH;^PT, including the effects 
of the finite lattice size [5]. The terms with coefficients Cj and c,- are analytic NLO countert- 
erms that cancel the renormalization-scale-dependence oi fi, and the terms with coefficients di^nnyp 
describe the leading effects of the non-zero lattice spacing. The functions /2,3 also depend on the 

S — T mass splitting A, which is included in the kinetic terms of Eq. (2.4). We set A = 200 MeV, 

(*) 
consistent with the l-l — At splitting from experiment and with our lattice data. 

To study the effect of the HYP smearing in the heavy-quark action on the scaling behavior, 

we performed initial fits that included all values of ?ihyp. and then successively removed the data 

with the largest values of nuyp- After excluding «hyp = 5, 10, the fits were stable and had good 

2-values. Our final results for the axial couplings are 

gi =0.449 ±0.047 Stat ± 0.019 syst, 

g2 = 0.84 ±0.20 Stat ±0.04syst, 

g3 = 0.71 ±0.12stat ±0.04syst. (4.5) 

The estimates of the systematic uncertainties in (4.5) include the effects of the following: NNLO 
terms in the fits to the a- and m^j-dependence (3.6%, 2.8%, 3.7% for g\, g2, g3, respectively), 
the above-physical value of the sea-strange-quark mass (1.5%), and higher excited states in /?;(?) 
(1.7%, 2.8%, 4.9%). The details of the analysis can be found in Ref. [6]. 

Figure 3 shows the pion-mass dependence of the fitted functions (g,)eft- The counterterm 
parameters c^ and c,- resulting from the fits are natural-sized (for /i = 4Kfji), and the NLO 
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contributions are significantly smaller than the LO contributions. We conclude that the SU{A\2) 
chiral expansion of (g,)eff convergences well for the pion masses used here. 

5. Summary 

We have calculated the heavy-hadron axial couplings using lattice QCD, including for the first 
time the baryonic couplings g2 and §3 in addition to the mesonic coupling gy. The analysis is 
based on data for the axial-current matrix elements at low pion masses, a large volume, and two 
different lattice spacings. We extracted gi,2,3 from this data by performing chiral fits with the full 
NLO expressions from HH;^PT [7]. As a consequence, the systematic uncertainties in our results 
(4.5) are much smaller than the statistical uncertainties. The numerical values of ^1,2.3 can be 
used to constrain chiral fits of lattice QCD data for a wide range of heavy-light meson and baryon 
observables. Furthermore, our results for the axial couphngs allow the direct calculation of certain 
observables in HH;t^PT, in particular the strong decay widths of heavy baryons [5,6]. 
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